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Abstract. The direct laser cooling of neutral diatomic molecules in molecular beams 
suggests that trapped molecular ions can also be laser cooled. The long storage time 
and spatial localization of trapped molecular ions provides the opportunity for multi- 
step cooling strategies, but also requires a careful consideration of rare molecular 
transitions. We briefly summarize the requirements that a diatomic molecule must 
meet for laser cooling, and we identify a few potential molecular ion candidates. 
We then perform a detailed computational study of the candidates BII+ and A1II+, 
including improved ab initio calculations of the electronic state potential energy 
surfaces and transition rates for rare dissociation events. Based on an analysis of 
population dynamics, we determine which transitions must be addressed for laser 
cooling and compare experimental schemes using continuous-wave and pulsed lasers. 
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1. Introduction 

The long-held notion that laser-cooling molecules is infeasible has been recently 
overturned by the transverse laser-cooling of SrF lU |2]. The direct Doppler cooling 
was possible due to the nearly diagonal Franck-Condon factors (FCFs) of the A-X 
transition and strong optical forces resulting from the short excited state radiative 
lifetime. Diagonal FCFs minimize the number of vibrational states populated by 
spontaneous emission during the cooling time. Similar to previous proposals [31 H], 
the number of relevant rotational states involved in the cooling cycle is minimized by 
judicious choice of the initial angular momentum state [1]. 

In ion traps, laser-cooled atomic ions can be used to sympathetically cool any 
other co-trapped atomic or molecular ion species O [6l [7j, and this technique is 
being used to study gas-phase atomic and molecular ions at very low temperatures 
in laboratories around the world. For example, sympathetically cooled ions have been 
used for mass spectrometry P [9], precision spectroscopy [IDl ttH [121 lEl ttH US], and 
reaction measurements [161 [13 [HI [El [201 [211 [22]. Because the impact parameter in 
ion-ion collisions is much larger than the molecular length scale, the internal states 
of sympathetically cooled molecular ions are undisturbed. A low-entropy internal 
state can be prepared either by producing the molecular ion through state-selective 
photoionization [23] or by taking advantage of the long lifetime of the sympathetically 
cooled molecular ions to optically pump the internal degrees of freedom [2H [251 [26] . 
Now that internal state control of sympathetically cooled molecular ions has been 
demonstrated, it is interesting to consider the possibility of eliminating the atomic 
coolant ions (and the accompanying laser equipment) and directly Doppler cooling 
certain molecular ions. Since molecular ions remain trapped for exceedingly long times 
independent of internal state, ion traps relax requirements on molecular transition 
moments and repump rates, thus offering a unique environment for molecular laser 
cooling. Direct cooling would also be desirable when the coolant ion may create 
unwanted complications (e.g., quantum information processing using trapped polar 
molecular ions coupled to external circuits [27]), or when the light driving the atomic 
cycling transition results in a change in the molecules' internal state (e.g., resonant 
absorption or photodissociation). Achieving the necessary cycling of a strong transition 
would also allow straightforward direct fluorescence imaging of trapped molecular ions 
at the single-ion level. 

On the basis of spectroscopic data available in the literature (see the appendix), 
our molecular ion survey found BH+ [28l [29l [30l [SH [32l [33] and A1H+ [HH [29l [35l [36l [37] 
to be among the most promising candidates. In this manuscript, we review all identified 
challenges of maintaining a closed excitation scheme for direct laser cooling of BH"*" and 
AlH"*" stored in ion traps. Unlike for the case of SrF, where slow vibrational decay relative 
to the interaction time allows a straightforward probabilistic approach to predict repump 
requirements [l], designing a laser cooling experiment for trapped ions requires careful 
modeling of vibrational decays within the ground state, since they result in diffusion of 
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parity and rotational quantum numbers. Additionally, for BH"*" and AlH"*" as with all 
molecular Doppler cooling experiments, rare decay and photofragmentation processes 
need to be considered, since they can occur on timescales comparable to Doppler cooling 
and much shorter than the ion trapping lifetime. 

This paper is organized as follows. In section |2] we discuss the molecular 
properties required for Doppler cooling. In section [3l we consider in detail two cooling 
candidates, BH"*" and AlH"*", present new quantum chemical calculations for these species 
(section IS.l.ip . discuss the calculation of the Einstein coefficients (section I3.1.2p . and 
present the rate-equation model used to determine the number of cooling photons 
scattered (section [3.1.3p . In section [321 we discuss the cooling scheme, including the 
effect of spin-rotation and spin-orbit splitting of the and ^11 states, respectively. 
The results of our simulation are in section \3.3\ and we present our calculations for 
rare events, which terminate the cycling transition, in section 13.3.31 In section HI we 
discuss the technological requirements to carry out the proposed experiments, and we 
propose the use of femtosecond lasers to provide multiple repump wavelengths from a 
single source. Finally, section [5] summarizes the concepts, simulations and prospects 
for Doppler cooling of molecular ions. In the appendix, we present a few additional 
classes of molecules that need to be studied in more detail to decide whether they are 
candidates for direct laser cooling experiments. 

2. General experimental considerations 

Doppler laser cooling requires the scattering of many photons, with each emitted photon 
on average carrying away an amount of energy proportional to the laser detuning from 
resonance [38] ■ To cool a two-level particle with a mass of ~ 10 amu and a visible 
transition linewidth of ~ 1 MHz from room temperature to millikelvin temperatures 
requires the scattering of 10^ — 10^ photons. Two-level systems are easily obtained in 
atoms, e.g., ^P-^S transitions in alkali metals (for neutrals) and alkaline earth metals 
(for ions). In both cases, the ground states have a closed shell with one valence electron, 
and the state of the atom is well described by the orbital and spin angular momentum 
of the valence electron [39] . 

Electronic two-level systems are impossible to obtain in molecules since the 
vibrational and rotational degrees of freedom introduce multiple decay paths not present 
in atoms. Spontaneous decay into non-cycling bound or dissociative states generally 
terminates the cycling transition long before the Doppler cooling limit is reached; 
however, a carefully chosen molecule can significantly reduce the probability of such 
decays. The ideal molecule would possess an excited and electronic ground state with 
equivalent potential surfaces separated by an energy offset in the optical region, resulting 
in perfectly diagonal FCFs. Most intuitive of the types of transitions which might 
lead to diagonal FCFs [ID], are transitions which excite a single electron from the 
highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital 
(LUMO). This transition will have negligible effect on the bond length when both the 
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HOMO and LUMO are either non-bonding or anti-bonding orbitals [H]. Alternatively, 
one may choose a molecule with an unpaired electron and a corresponding hole in one of 
the highest bonding orbitals, in which the optical transition moves the hole to another 
bonding orbital. 

In addition to considering the FCFs, one must take into account the effect of 
rotational branching due to spontaneous emission, unexpected decay paths due to the 
break down of the Born-Oppenheimer approximation and dipole-forbidden transitions. 

2.1. Diagonal Franck- Condon factors 

When a decay from the excited state to the ground state occurs, the branching ratio into 
final vibrational states is determined from the Franck-Condon overlap, which measures 
the similarity between the excited and ground vibrational wavefunctions [42j. Molecules 
which have similar excited-state and ground-state potential energy curves (PECs) reduce 
the total number of repumping lasers required, so it is useful to determine the repumping 
requirements as the two PECs are made more dissimilar. 

High-accuracy quantum-chemical calculations of the excited PEC are not typically 
available; however, the Morse potential, with only three parameters, is a reasonable 
approximation for many bound electronic states. This approximation is particularly 
true near the bottom of a potential well, which is energetically distant from the region 
of non-adiabatic couplings that alter the dissociative asymptotes of the potentials. 
Assuming excited vibrational levels have infinite lifetimes, the number of transitions 
that are required to approximately close a cycling transition depends strongly on all 
of the three parameters. Figure [1] shows the effect of changing the equilibrium bond 
distance and vibrational frequency for fixed anharmonicity. The number of vibrational 
levels that must be addressed to have less than 1 part-per-million population loss is very 
sensitive to differences in the ground and excited state bond lengths. In the ideal case 
of perfectly diagonal FCFs, a single laser is required; a difference of only 5% in the bond 
lengths requires from four to seven lasers. 

2.2. Elimination of rotational branching 

In the dipole approximation, each emitted photon can change the total angular 
momentum of the molecule by 0, or ±1 h. The emission results in the population of 
additional states through Stokes and anti-Stokes processes. Each additionally populated 
angular momentum state will need to be addressed in order to achieve a high scattering 
rate and efficient laser cooling. 

The number of angular momentum states populated depends on which Hund's case 
(a) or (b) is important |l2]. For the purpose of our discussion here, we will initially 
ignore spin (equivalent to Hund's case (b)) and later look at complications which occur 
in Hund's case (a). We will also ignore nuclear spin and the resulting additional hyperfine 
states and transitions (the resultant hyperfine splitting is sufficiently small that in 
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Figure 1. The effect of excited state vibrational constant and equilibrium bond length 
on the number of optically pumped vibrational states required to have less than 1 
part-per-million population loss after spontaneous emission. Calculations are based 
on probabilistic analysis of FCFs of an hypothetical A ~ X system. White space 
represents parameter regions in where at least 7 vibrational bands must be addressed. 
The ground and excited state PECs are represented by Morse functions with a fixed 
anharmonicity of WgXe of 50 cm~^. The ground state, X, has an equilibrium distance 
Re of 1.32 A and a vibrational constant of 2750 cm"-'^. Laser cooling with one laser 
is possible only in a narrow region (in black) where the ^-state spectroscopic constants 
are similar enough to i?e and Wg. 



practice these additional states can easily be addressed by use of an electro-optical 
modulator (E0M)[I1 [2]). 

A— doubling breaks rotational symmetry about the bond-axis, resulting in electronic 
states of well-defined symmetry relative to the plane defined by the axis of nuclear 
rotation and the bond-axis. A key feature of the SrF experiment [H [2] is pumping from 
the ground state f'=0, N'=l, K'=l state with negative parity to the excited f=0, N=0, 
K=l state with positive parity, where N is the rotational angular momentum of the 
molecule, and K is the total angular momentum without spin. As a result, decay from 
the excited state is limited to those v allowed by Franck-Condon factors and only via 
the Q branch. The desired transitions and the expected spontaneous emission channels 
are shown in figure El 

This method requires the control of the initial angular momentum population, as 
not all the molecules in the experiment start out with the desired or J. However, 
the preparation of an specific initial state of the molecular ion can be achieved by laser- 
cooling the internal states [211 |25l [26] , by cooling the internal states through interaction 
with cold neutral atoms |13], using a cryogenic ion trap [S] to limit the number of 
energetically accessible states and the rate at which they mix, or using state-selective 
photoionization of a neutral molecule |23] . 
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Figure 2. Energy level schematic of proposed ^II-^E transition. Red dashed levels 
are odd parity and blue solid levels are even parity. K refers to the total angular 
momentum without spin, and v and v' label the vibrational levels of the excited and 
ground states. Starting in the K — 1 state and exciting the Q-branch (solid arrows) 
results in a closed transition modulo the vibrational states (wiggly arrows). 

Additional rotational states can be populated by blackbody redistribution between 
different rotational states, or by the decay of vibrationally excited states through 
spontaneous emission. This process does not occur for homonuclear diatomic ions. 
For heteronuclear diatomic ions, blackbody redistribution results in an increase in the 
number of required repumps if the redistribution occurs on cooling timescales. For BH"*" 
and A1H+ we find that rotational redistribution requires tens of seconds, which is longer 
than the cooling time, and is included in our simulation (see section [373]) . Vibrational 
radiative relaxation can limit the number of vibrational bands that must be addressed, 
at the cost of repumping from a larger set of rotational states due to the diffusion in J 
during the vibrational cascade. 

If vibrational relaxation is sufficiently fast, a probabilistic FCF approach to 
determine repumping requirements is no longer valid, since the FCF approach assumes 
vibrational levels are infinitely long-lived. For example, if we cycle on the f ' = -H- f = 
transition, and repump on the v' = 1 ^ v = transition, the ion will occasionally decay 
via the f = — )■ = 2 — )■ f ' = 1 and the f =' 1 — )■ f ' = channels, resulting in parity 
flips and occupation of higher rotational states. 

A simple figure of merit (FOM) for the probabilistic FCF approach is calculated by 
comparing the total time required to Doppler cool the molecule to the vibrational decay 
rate, since all levels are approximately equally populated at steady-state. The FOM 
is given by j^^, in which Aqqi is the excited-state decay rate, is the total number 
of photons scattered to Doppler cool, and Aiiqi is the vibrational decay rate. Using 
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an estimate of ^ 10^ we calculate a FOM of 1, 20, 250 for BH+, A1H+, and SrF, 
respectively. We expect that the effect of vibrational decay is most significant for BH+, 
followed by AlH"*", which have much smaller FOMs than does SrF. Our rate-equation 
simulations (see section 13. 3p confirm this prediction, with vibrational decay in BH"*" and 
AlH+becoming important at short and long times, respectively, and vibrational decay 
in SrF unimportant on cooling timescales. 

2.3. Other decay paths 

The discussion so far has relied on the Born-Oppenheimer approximation and only 
considered dipole-allowed electronic transitions. These approximations reliably describe 
the molecules of interest for timescales up to a few microseconds. However, long ion-trap 
lifetimes require the consideration of slow processes driven by violations of the model. 

Ideally, the A state is stable against dissociation. However, constant excitation of 
the molecule allows for rare events to occur, such as dissociation via the repulsive part of 
the ground state potential. This dissociation can be caused by spin-orbit coupling or L- 
uncoupling, which violates the Born-Oppenheimer approximation by mixing the ground 
and excited electronic state potentials |15]. In the case of L- uncoupling, it is possible to 
pick excited states such that symmetry effects forbid dissociation from occurring. Spin- 
orbit coupling also exists, but the calculated dissociation rate ranges from kHz to Hz 
and, to our knowledge, has never been measured. (In most molecular beam experiments 
these rare events are unobserved, since the dissociation timescales are long compared 
to laser- interaction timescales.) It is also possible that cooling or repumping lasers can 
cause photodissociation by coupling to the A state to a higher-lying dissociative state. 

Additionally, transitions that are not electric dipole allowed can occur due to mixing 
of parity states by electric fields. For a compensated ion trap, the stray DC electric 
field at the axis can be effectively cancelled |16]. Population of opposite parity states is 
inevitable in the long-time limit, due to magnetic dipole transitions (those of A J = 0, ±1 
terminate on states of the wrong parity). These transitions occur at a rate of 1 transition 
or less in 10^, also requiring attention when forming a highly-closed cooling cycle [3]. 

Figure [3] graphically summarizes all the mechanisms by which the diatomic 
molecular ion can leave the cycling transition, including optical photon emission to an 
excited vibrational state followed by the emission of an infrared photon that transfers 
population to a lower vibrational state. As mentioned in the previous subsection, while 
this process leads to diffusion in the rotational states and introduces states of additional 
parity, it may ease the experimental requirements by reducing the number of vibrational 
bands that must be addressed. 

3. Example laser cooling cycle for BH+ and A1H+ 

The A^H-X^E+ system of both BH+ and AlH"*" are well studied in emission from hollow 
cathode discharges SHI EO], and chemiluminescent ion-molecule reactions ^91 [35] 



Challenges of laser-cooling molecular ions 



8 




Figure 3. Energy level schematic of proposed ^II-^S transition showing transitions 
that result in either dissociation or occupation of a rotational state outside of the closed 
transition. 

have yielded detailed rotational constants for the first two vibrational levels of the X^S"*" 
ground state along with low-resolution spectra of a few diagonal vibrational bands. For 
the BH^ cation, 15 rovibrational energy levels of the ground state (v=0-4 and N=0- 
2) have been assigned by the extrapolation of photo-selected high Rydberg series of 
neutral BH [32J. During those experiments, a very strong optical excitation of the ion 
core was observed when the energy of the scanning laser matched the A-X transition of 
the cation |50]. Photon absorption rates observed during this isolated core excitation 
were consistent with Einstein transition probability coefficients calculated theoretically 
by Klein et al. [29]. Another favorable property of this molecular ion is that the first 
quartet states are energetically above the A-X system [311 [51] , eliminating the possibility 
of intervening electronic states of different spin to which the upper state could radiate 
and terminate the cycling transition. With all this information, it is possible to consider 
a laser cooling experiment for trapped BH"^ ions. Similarly, we expect A1H+ to share 
most of these favorable characteristics. 

3.1. Methods 

3.1.1. Ah initio calculations for BH^ and AlH^ To study the processes outlined in 
section [21 potential energy functions and dipole moments for the X^S'^, A^II, B^S"*", and 
a^n states of BH+ and A1H+ are required. Additionally, the transition dipole moments 
and spin-orbit coupling matrix elements between these states are needed. For BH"*", it 
is possible to solve the electronic Schrodinger equation exactly (within a basis) through 
full configuration interaction (FCI). By constraining the Is orbital of B to be doubly 
occupied, a large aug-cc-pV5Z basis [52] can be used to obtain the potential energy 
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functions. The spectroscopic constants obtained at this level of theory [FCI(3e~)/aug- 
cc-pV5Z] accurately reproduce the experimentally determined values (shown in tabled]), 
and were computed by fitting a fourth-order polynomial to five energy points centered at 
Re and evenly spaced by 0.005 A. For AlH"*", FCI was not applicable because correlating 
the important 2s and 2p electrons of Al leads to too many determinants in the expansion 
of the wavefunction. Coupled-cluster methods and their equation-of-motion variants 
provide a practical alternative |53l [51]. In particular, CCS and E0M-CC3 [55] were 
applied with an aug-cc-pCVQZ basis set [52]. Again, the Is orbital was constrained to 
be doubly occupied. The FCI and E0M-CC3 computations were performed with PSI 
3.4 [56]. 

Dipole moments, transition dipole moments, and spin-orbit coupling matrix 
elements were obtained from multireference CI wavefunctions (MRCI). The MRCI 
wavefunctions add single and double excitations to a state-averaged complete active 
space self-consistent field (SA-CASSCF) reference [SUES]. The SA-CASSCF orbitals 
were optimized with an active space consisting of three electrons in the 2al7Tx^'n'y3a4a5a 
and 4cr27r^27rj^5cr6cr7cr orbitals for BH+ and AlH"*", respectively. In the MRCI 
wavefunctions, single and double excitations from the Is B orbital and from the 2s 
and 2p Al orbitals were allowed. Again, the aug-cc-pV5Z and aug-cc-pCVQZ basis sets 
were used for BH"*" and A1II+, respectively. The MRCI computations were performed 
with MOLPRO [59]. 

All of the PECs, dipole moments, transition dipole moments and spin-orbit 
matrix elements (doublet-doublet, quartet-doublet and quartet-quartet) generated by 
the calculations described in this section have been included as a supplementary file. 

3.1.2. Einstein A and B coefficients The Einstein A coefficients connecting two states 
are calculated using 

A _ ^u,l 9u,l Q ^2 /i\ 

- 3^2XTT'^^''^"'^' 

where uJu^i/2t: is the transition frequency, c is the speed of light, eo is the permittivity 
of free space, gu,i is a degeneracy factor, Ju is the total angular momentum of the upper 
state, Su,i is the Honl-London factor, and Vu^i is the transition dipole-moment matrix 
element. The B coefficient is calculated from the A coefficient using 

B^, = (2) 

3.1.3. Rate-equation model A rate-equation approach was used to model the 
population dynamics of BH"*" and AlH'*' allowing determination of the average number 
of Doppler-cooling photons scattered before a molecule was pumped into a dark state 
(ignoring predissociation and photodissociation described in section r3.3.3p . A similar 
approach has been used for internal-state cooling [24l[26l|25]. In particular, we solved 

^^MP, (3) 
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Table 1. Spectroscopic constants for the X^I]+ and A^II states of BH+ and A1H+. 
Comparisons to previous theoretical and the most precise experimental values, where 
available, are also given. Re is the equilibrium bond length, B^., D^., and ae 
parameterize the rotational energy level spacings [42j, We and ujeXe are the harmonic 
and anharmonic frequencies of the potential at equilibrium, and Dq is the dissociation 
energy. Units are cm~^ unless otherwise specified. 
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where P is a vector consisting of the rovibrational levels we include in our model, 
ordered by increasing energy, and M is a x vector consisting of the Einstein A and 
B coefficients which couple the different rovibrational states. Explicitly, the population 
of a given state follows 

^p, j=i-T- i=«-l j=N 

j=l j=l j=i+l 

j=N j=i-i j=N 

+ E ^J^P3 + E ^M^3^)P3 + E ^M^,i)Pr (4) 
j=i+l j=l J=«+l 

Here, the i^^ and j*'^ states are connected by Einstein coefficients denoted by Aij, 
Bij, and Bji which correspond to spontaneous emission, stimulated emission, and 
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Figure 4. Born-Oppenheimer potential energy functions of the first two doublet 
electronic states of BH+ at the FCl(3e~)/aug-cc-pV5Z level. Black dots are the actual 
FCI calculation results, and the solid lines are splines with analytical functions as 
described in the LEVEL 8.0 and BCONT 2.2 manuals [Ml ES]. Exact numerical 
solution of the radial Schrodinger equation yields the plotted wavefunctions, which are 
used to calculate the various decay and emission rates described in the text. 

absorption, respectively, and p{uJij) is the spectral energy density at a given frequency, 
ojij. The Einstein coefficients were calculated using the potential energy curves and 
dipole moments discussed in section I3.1.1I 

The average number of Doppler cooling scattered photons was calculated by 
numerically solving equation (j3]), and multiplying the population in the excited state 
of the cycling transition by its spontaneous emission rate. Counting photons emitted 
into the ground state from the excited state, regardless of whether the excited state 
was populated by a cycling-laser photon or a repumping-laser photon, is equivalent 
to counting the number of cooling-cycle absorption events. A simple estimate of the 
number of scattering events required for cooling can be obtained by considering the 
average energy lost per scattering event. For a laser detuning r/2 from resonance, where 
r is the natural linewidth, then the average energy lost per scatter is /S.E = hV /2. The 
A2n-X2E+ linewidths of BH+ and A1H+ are T = 27r x 0.7 MHz and T = 27r x 2.6 MHz, 
respectively. Starting from T = 300 K, cooling to mK requires ~ 10^ scattering 
events for both BH"*" and AlH"*" for our laser detuning. We may also take advantage of 
the Doppler width at higher temperatures by detuning the laser further from resonance 
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Figure 5. Energy level schematic showing cooling and repumping transitions for the 
case of negligible vibrational and rotational relaxation within the excited and ground 
state manifolds (not to scale). Ground and excited states are labeled \g^, j,) and j), 
respectively. 

to increase the energy lost per scattering event, reducing the laser detuning as the ions 
are cooled. Using this method, we estimate a lower limit of ~ 10^ scattering events 
are required to cool to mK. 

Finally, it should be noted that a rate-equation approach does not account for 
coherence effects which can result in dark states. For example, in atomic systems with 
a low-lying D state, two lasers are typically used for Doppler cooling to address an 
5 — 7- P (cooling) transition and a D — t- P (repumping) transition. When both lasers 
are detuned from resonance by similar amounts, dark resonances occur and cycling 
terminates [66]. The cooling and repumping transitions in our molecular system also 
exhibit a similar A-system (see figures [5] and [6]) . Experimentally, these dark resonances 
can be avoided either by explicitly varying the laser detuning, or by shifting the energy 
levels by different amounts with an externally applied time- varying field [661 167] . 

3.2. Cooling scheme 

The general cooling scheme, ignoring spin, was discussed in section [2] and figure El 
The inclusion of spin results in spin-rotation splitting for the states and spin-orbit 
splitting of the ^11 states [12]. We label the X^E+ and A^ni/2 states as \g^i'j') and |e^j). 
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Figure 6. Energy level schematic showing cooling and repumping transitions for the 
case of non-negligible vibrational and rotational transitions within the excited and 
ground state manifolds (not to scale). 

respectively. The Zeeman substructure is explicitly ignored in our model; however, the 
resulting reduction in the overall scattering rate due to differing multiplicities of the 
ground and excited states is included as a degeneracy factor in the Einstein A and B 
coefficients. In practice, an externally applied magnetic field of a few Gauss is sufficient 
to lift the degeneracy, and the laser polarization can be chosen to ensure that all Zeeman 
sublevels are addressed [66l |67]. Additionally, both AlH"*" and BH+ have hyperfine 
structure, which is ignored in our model. 

We assume that population begins in the \gQi/2) state and a laser drives transitions 
to the Icq ^^2) state. From the excited state, population may decay into either the \g~, ^^2) 
state or the \g~, 3^2) state, with decreasing probability for increasing v'. The cooling 
transition, corresponding to v' = 0, requires driving both the \gQi/2) < — > kdi/2) 
l5'o'3/2) ^ — ^ 1/2) transitions, as shown in figure |5] with labels CI and C2. The ground 
states belong to the same rotational level {K' = 1), so this cooling scheme is equivalent 
to the one discussed in section 12.21 

If vibrational levels with f ' > 1 are long-lived, then decay between vibrational levels, 
which result in parity flips, can be ignored. In this case, we drive \g~, ^^2) ^ — ^ kj-i 1/2) 
and \g~/ ^/2) < — ^ I^J-i 1/2) transitions to pump population back into the cycling 
transition, as shown by the Rl, R2 and unlabeled arrows in figure [H 

However, when decay between vibrational levels is fast, the resulting parity flips 
may populate the \g^f^ji) states, so lasers connecting these states to an excited state 
must also be included to avoid population buildup in the even-parity states. Multi-stage 
vibrational decay also results in angular-momentum diffusion. For instance, decay from 
the [cq ^^2) populates the \gi^/2) state, which decays into the \gQ^/2) state. Population is 
pumped out of states with J' > 3/2 by driving \g^j/) < — > \eQj,_i). This cooling scheme 
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Figure 7. Energy level schematic showing rotationless lifetimes of different BH+ 
levels (not to scale). Lifetimes for A-^>X transitions are denoted by r and lifetimes for 
transitions within the ground state are denoted by r'. All lifetimes are calculated using 
the ab initio potential energy curves and dipole moments discussed in section [3.11 

is illustrated in figure El 

In figures [7] and [9] we show the rotationless lifetimes for BH+ and A1H+ respectively. 
Since both AlH"*" and BH"*" have relatively fast decay rates (up to 100 s~^) between 
vibrational levels, we find that they require cycling and repumping on both parities. 
In the following section, we examine the details of each molecule separately, and 
determine their repumping requirements. Predissociation and photodissociation rates 
are considered in the last subsection. 

3.3. Simulation results 

3.3.1. BH^ In figure [8] the total number of scattered photons for different laser 
configurations is plotted. Population is initially in the \gQii2) state. The simulation 
includes dipole-allowed transitions between electronic states, and vibrational and 
rotational transitions within the X^E"*" state. Blackbody radiation is also included in 
the model, for a temperature of 300 K. In Case (a-i), two wavelengths are used to 
drive the odd-parity cycling transitions: CI and C2 in figures [5] and |6l Population is 
quickly pumped into the v' = 1 manifold, and a total of ~ 70 photons are scattered. 
In Case (a-ii), an additional two wavelengths are introduced to drive the even-parity 
cycling transitions: C3 and C4 (see figure |6]). Again, population is quickly pumped into 
the v' = 1 manifold; however, at longer times (t ~ t^q), population leaks back into the 
even-parity ground states and is driven by the even-parity cycling transitions, resulting 
in a total of A^ ~ 380 scattered photons. The delay between successive photon count 
plateaus is due to population building up in the v' = 1 manifold. 
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Figure 8. Plot of the integrated cycling-transition photon count for the BH+ 
simulation using different laser configurations, driving subsets of transitions shown 
in figure [51 The hatched region corresponds to upper- and lower-limits for Doppler 
cooling to mK temperatures, corresponding to fixed and swept detuning, respectively. 



Case (b-i) includes the transitions from Case (a-i) plus an additional two 
wavelengths to drive the odd-parity repumping transitions: Rl and R2 in figures [5] 
and [6l ~ 8000 photons are scattered, with the dominant dark-state leak being 
the even-parity ground states. Case (b-ii) extends Case (a-ii) with an additional four 
wavelengths to drive even- and odd-parity repumping transitions: R1-R4. A total of 
^ 21000 photons are scattered, which plateaus at just above the swept- detuning 
threshold. 

As predicted in section [3^ vibrational decay plays a significant role in depopulating 
the higher-lying vibrational ground states, resulting in parity fiips and rotational 
diffusion. In Case (b-ii), the parity-fiip issue is addressed, so population built up in 
the 1(7^5/2) states. Case (c-ii) includes the same wavelengths as Case (b-ii) and also two 
additional wavelengths to drive odd- and even-parity P-branch repumps: PRl, PR2 (see 
figure E]). The increase in counts (A^ ~ 600, 000), which is well above the swept-detuning 
threshold required for significant cooling from room temperature. 

3.3.2. AlH^ The simulation of A1H+ is plotted in figure [TOl Again, we assume that 
population starts in the |5'(7i/2) state. The simulations include transitions between 
electronic states, within the ground electronic state, and include a blackbody source 
of radiation at 300 K. In Case (a) four wavelengths are used to drive only the cycling 
transition: C1-C4. A total of A^ ~ 200 photons are scattered. 

In Case (b) an additional four wavelengths are used to drive the repumping 
transitions: R1-R4. There is a significant increase in counts (A^ ~ 15000), and a 
noticeable deviation from the results of BH+. Ai t ^ IQ'^ s, there is a reduction in 
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Figure 9. Energy level schematic showing rotationless lifetimes of different A1H+ 
levels (not to scale). Lifetimes for X^A transitions are denoted by r, and lifetimes 
for transitions within the ground state are denoted by r'. All lifetimes are calculated 
using the ab initio potential energy curves and dipole moments discussed in section 

EH 




10"'' 10'^ 10''' 10'^ 10"^ 10"' lo" lo' 

time (s) 



Figure 10. Plot of the integrated photon count for the A1H+ simulation using different 
laser configurations, driving subsets of transitions shown in figure [6] The hatched 
region corresponds to upper- and lower-limits for Doppler cooling to mK temperatures, 
corresponding to fixed and swept detuning, respectively. 
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rate, resulting from population building up in the v' = 2 level. For BH"*" the rate of 
populating f ' = 2 is 55 times slower than the decay from f ' = 2 to f ' = 1 {tq2/'T2i = 55). 
For AlH"*" the rate of populating f ' = 2 is 500 times faster than the decay from v' = 2 
to f ' = 1 (T02/T21 = 1/500), resulting in population trapping. 

In Case (c) an additional four wavelengths are applied to drive the repumping 
transitions: R5-R8. This results in a total of ~ 210, 000 photon counts, which is 
well above the swept-detuning cooling threshold. Interestingly, the second vibrational 
repump is sufficient for cooling of A1H+ whereas BH"*" requires optical pumping of the 
{95/2) states. 

In Case (d) AJ = —1 repumping transitions are included: PRl and PR2. The 
resultant photon count {N ^ 700,000), which is similar to the number of counts we 
obtain in Case (c-ii) of BH"'' . We find that the results of AlH"^ are similar to those 
of BH+ in that they each require repumping of both parities, due to vibrational decay; 
however, AlH"*" differs in that it requires an additional vibrational repump, while not 
requiring additional rotational repumping. 

3.3.3. Calculation of predissociation and photofragmentation rates The cooling scheme 
requires a continuous repopulation of the vibrational ground level of the A^II excited 
electronic state. One possible mechanism that stops the cooling cycle is isoenergetic 
predissociation of the molecular ion by the coupling of the bound excited state to a 
dissociative state of the X^E"*" repulsive wall. Bound and continuum wave functions can 
interact via spin- and rotational-orbit couplings. The expression for such predissociation 
rates, in units of s~^, is obtained from Fermi's golden rule [61] 

h{v,J) = '^\{4je,j'{R)\Ms{R)\^.,j{R))\\ (5) 

where Ms{R) is a term in the molecular Hamiltonian, which is neglected in the Born- 
Oppenheimer approximation and describes the coupling between the initial and final 
electronic states. The lipyjlR)) initial-state wavefunction is space-normalized and the 
\ipE,j'{R)) wave function of the continuum is energy-normalized [IS]. The calculation 
uses a one- dimensional density of states in the bond coordinate to determine the 
energy-normalization, d{E) oc ^ fi/{E — D), where is the reduced mass and D is 
the dissociation energy of the final electronic state |64j . 

The predissociation depends strongly on the slope of the repulsive part of the 
dissociative potential. In order to estimate the uncertainty in the predissociation 
calculations, we first calculate the overlap integral of the wavefunctions assuming 
and multiplying by a bond-length independent spin-orbit coupling of 13.9 cm"^ [3U] . 
For BH+, by varying the details of the bound-state wavefunction, we find that the 
predissociation rate varies by more than four orders of magnitude from 10^ to 10~^ 
s~^. Representing the A^n(v=0) wave function by a symmetric Gaussian from the 
solution of the harmonic oscillator, the rate results in the most disastrous scenario, 
4.4 X 10^ s~^. Using the solution of a Morse potential gives 6.2 x 10~^ s~^. Combining 
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the continuum and bound wavefunctions computed using LEVEL 8.0 and BCONT 2.2, 
respectively, and by inputting the FCI(3e~)/aug-cc-pV5Z A^H and X^S+ PECs, the 
predissociation rate is 1.6 x 10~^ s~^ (potentials and wavefunctions shown in figure H]). 
Including the spin-orbit coupling function calculated in section I3.1.H as prescribed in 
equation (jS]), the predissociation rate results in a value of 0.03±0.02 s~^. For AlH"*" 
(using E0M-CC3 potential energy functions and MRCI spin-orbit coupling matrix 
elements), the predissociation rate is 0.2±0.1 s~^. The quoted uncertainties come from 
interpolation noise in both the repulsive part of the X^E+ state and the point piecewise 
polynomial used to fit the R-parametrized spin-orbit couplings. This uncertainty due 
to interpolation is negligible for the calculation of optical decay rates (~10 parts per 
million) reported elsewhere in this manuscript. 

Another problematic event is the excitation from the A^n state to either bound or 
dissociative B^E+ levels, by absorption of a cooling or repump photon. For BH"*", PECs 
from two different ab initio levels of theory place the dissociation limit, Vum, either 
below or above the energy accessible by the 379 and 417 nm photons exciting from the 
A state (see table [2] and figure [TTj) . Although both results agree with predictions based 
on experimental observations (due to a 2000 cm~^ uncertainty), these calculations yield 
two different photon absorption scenarios. The first (MC-SCF) would photodissociate 
the molecule, and the second (FCI) would steal population from the cooling cycle to a 
bound vibrationally excited level of the B2S+ state. For A1H+, MC-SCF and E0M-CC3 
predict Vum to be below the total energy accessed by any of the three required cooling 
and repump wavelengths (360, 381 and 376 nm). 

The photodissociation cross section (in cm^ /molecule), is [68] : 

(6) 

S is the usual Honl-London rotational intensity factor (Sj' /{2J + 1) is set to 1 in this 
section), and Mg{R) is the transition moment function, including the ratio of initial to 
final-state electronic degeneracy factors. Allowed values of J' are given by the usual 
rotational selection rule. Similar to equation ([5]), this expression assumes that the 
continuum radial wavefunction amplitude is energy normalized. 

Assuming saturation intensities, calculated from Einstein A coefficients for the 
main cycling and repump transitions, we calculate the photon fiux for each transition. 
Is (T), listed in table [2j The photodissociation rates are the multiplication of these 
laser intensities by the cross sections calculated from equation ([6]). Continuum and 
bound wavefunctions were calculated using BCONT 2.2 and LEVEL 8.0, respectively, 
and using the E0M-CC3 potential energy curves and the MRCI transition moment 
functions. We perform simulations using the MC-SCF potential energies and transition 
dipole moments from [29] for comparison. We determine that photodissociation occurs 
on a timescale of tens of hours for BH"*", and several days for AlH"*". 
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Figure 11. Sensitivity of problematic two-photon channels to precise details of 
molecular structure. Potential energy functions of the A = 1 electronic states 
of BH+ are plotted at the a) MC-SCF level from [2g, and b) FCI(3e-)/aug-cc- 
pV5Z. Black dots are the actual ab initio results and the solid lines are splines 
with analytical functions described in the BCONT 2.2 and LEVEL 8.0 manuals [Ml 
165] . The wavefunctions shown, which are used to calculate photofragmentation 
or photoabsorption rates, are positioned at the eigenvalues that solve the radial 
Schrodinger equation. The length of the arrows represents the actual experimental 
energy of the CI and Rl photons. 

For BH"*", the FCI potential energy curves predict that the cycling and repump 
wavelengths result in an excitation from the A^II state into a bound level of the B^S"*" 
state. We believe that the probability of resonant absorption to the B^S"*" state is small 
due to the narrow laser linewidths relative to the vibrational energy spacing and due to 
the poor Franck- Condon overlap. The most intense B — A band reaching this energy 
region of the B^S+ state is the (11, 0) with the value for A of 1.6 x 10^ s~^. Population 
of this state results in vibrational diffusion which is not repumped in our cycling scheme. 

Given the 2000 cm~^ experimental uncertainty in the dissociation energy, we do 
not have enough information to determine if the photoabsorption (predicted by FCI) or 
photodissociation (according to MC-SCF) processes will present a problem for Doppler 
cooling 

4. Prospects for Cooling 

Our simulations show that we require a total of eight wavelengths to Doppler cool BH"*" 
and AlH+(see figures [H]and[l0]); however, cooling does not necessarily require eight 
unique lasers. For both ions, two of the cycling transitions, CI and C2, can be addressed 
by a single laser, since these two wavelengths differ only by the spin-rotation splitting 
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Table 2. Energies and rate coefficients relevant to possible pliotofragmentation events. 
huj (T), A (T), and Is (T) correspond to the photon energy, A Einstein coefficient, 
and saturation intensity for each A^n-X^S"*" transition T as labeled in figures [5] and 
[S] The BH+ cooling schemes do not involve the transition R5. M"''('^P)-M+(^S) is the 
energy difference between the lowest singlet and triplet energy states of B+ or A1+. 
This should correspond to the energy difference between the dissociative asymptotes of 
the X^S+ ground state and the common A^II and B^E+ excited states Yum- Etotai is 
the energy of the A^II (ti = 0, J = 1/2) state plus the photon energy of the transition 
T. Photofragmentation does not occur in the FCI(3e~)/aug-cc-pV5Z) calculation for 
BH+ (see text), but the X^S+ De does not preclude either possibility due to current 
experimental uncertainties. 
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(Ausr ^ 540 MHz for BH+ ^ and Ausr ^1.7 GHz for A1H+ This splitting 

can be achieved from a single laser source by using either an acousto-optic modulator 
(AOM), or an electro-optic modulator (EOM). The other two cycling transitions, C3 
and C4 each require their own source since their splitting is large {Au ~ 1 THz for BH+ 
and Au ^ 2 THz for A1H+). The cycling transitions require CW lasers that are narrow 
relative to the electronic linewidth. 

BH"*" requires driving the transitions Rl-4, corresponding to repumping v' = 1. 
These transitions also span 1 THz, so they would require a total of three CW lasers, 
since spin-rotation splitting between R2 and R3 can be bridged by an AOM. However, 
repumping on all four transitions can be achieved with a single femtosecond pulsed 
laser, with a typical bandwidth of 100 cm~^(3 THz). Using femtosecond repumps 
generally requires optical pulse shaping, in order to avoid R-branch (rotational heating) 
transitions. Pulse-shaped pumping has been used to achieve vibrational cooling of Cs2 
molecules, with a resolution of 1 cm~^ (30 GHz) [69] - which is sufficient resolution 
for rotationally-resolved repumping of hydrides. A femtosecond laser can also be used 
to drive the P-branch repumping transitions PRl and PR2; however, these transitions 
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would require a second laser since the Rl-4 transitions differ from the PRl and PR2 by 
more than 100 cm~^ (3 THz). It should noted that repumping of PRl and PR2 may 
not be necessary, since the total photon count exceeds the number required for Doppler 
cooling in the swept-detuning case (see figure [8]). 

AlH"*" requires driving transitions Rl-4, corresponding to repumping v' = 1 and R5- 
8, corresponding to repumping v' = 2. The repump Rl-4 span 2 THz, which is within 
the bandwidth of a femtosecond laser. Transitions R5-8 also span 2 THz; however they 
do not overlap with transitions Rl-4 {Au ~ 9 THz), so they would require a second 
femtosecond laser. 

If a broadband source for repumping is used in conjunction with CW lasers for 
cooling, the laser requirements for BH+ and A1H+ are reduced to five and six lasers 
driving 10 and 14 different transitions, respectively. For BH+, repumping transitions 
R1-R4 require one femtosecond laser, and repumping transitions PRl-2 require a second 
femtosecond laser. For A1H+, repumping transitions R1-R4 require one femtosecond 
laser, and repumping transitions R5-R8 require a second femtosecond laser. The PRl-2 
wavelengths would require a third femtosecond laser but they are not strictly necessary 
(see figure fTOj) . 

The comb-like nature of the frequency components of a femtosecond laser can be 
problematic if a transition which we wish to drive happens to fall between two comb 
lines (typically separated by 80 MHz). Various modulation techniques, either internal 
or external to the femtosecond laser, can be used to address this complication. 



5. Conclusions 

The key element for direct laser cooling of molecules [HI HI [H [2] and molecular ions 



are diagonal FCFs. From our survey of a number of molecules (see [Appendix AD 
we determined BH+ and A1H+ to be favorable candidates since a probabilistic FCF 
analysis (excluding vibrational decay) suggested that three and two vibrational bands, 
respectively, needed to be addressed to scatter 10"^ photons. A detailed analysis of 
state population dynamics (including vibrational decay) has shown that BH"*" requires 
addressing two vibrational bands while AlH"*" requires addressing three. For both 
molecular ions we find that vibrational decays result in diffusion of both parity and 
angular momentum; thus, implementation of this cooling scheme requires addressing 
the molecules with a large number of repump wavelengths. 

To minimize the total number of lasers required, a femtosecond laser can be used, 
since multiple transitions fall within the laser bandwidth. However, for many of these 
transitions, the frequency of excitation needs to first be experimentally determined with 
higher accuracy. The spectroscopy can be performed and the direct cooling technique 
optimized by initially trapping the molecular ions with laser-cooled atomic ions and 
monitoring the sympathetic heating and cooling of the atomic ion by the molecular ion 
|70] . Ultimately, dissociation of the molecular ion represents a fundamental limit for 
continuous Doppler cooling of the species presented here. 
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The attempted laser-cooling of BH"*" and AlH"*" will either result in the first directly 
laser-cooled molecular ions and/or and accurate measurement of rare predissociation 
events. Either way it will improve our knowledge of the spectroscopic splittings of BH"*" 
and A1H+ by an order of magnitude. Together with experimental dissociation rates, 
measurements of these splittings yield the most stringent test of advanced electronic 
structure theories; as the five-valence-electron open-shell diatomic BH"*" is among the 
simplest stable molecules and serves as a prototype for the development and testing of 
such theories [63l|28l[29l|3lll7Il[72l|6llini[T3l[62]. 

More generally, any molecular ion can be sympathetically cooled, and the long 
trapping hfetime allows slow processes to be investigated. Thus, sympathetically cooled 
molecular ions could serve as a valuable tool in determining the rates of higher-order 
processes such as predissociation and photofragmentation, which are of concern for direct 
Doppler cooling applications. 
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Appendix A. Potential Candidates 

Doppler cooling candidates are divided into two main categories. The first category 
corresponds to transitions in which the electron is excited from the HOMO to the LUMO, 
and the second category corresponds to transitions in which a hole moves from one 
orbital to another. Each category is subdivided according to grouping on the periodic 
table. Groups which we discuss are chosen, based on chemical intuition, to form open- 
shells with vertical excitations that do not perturb the chemical bond. Probabilistic 
predictions in this appendix are based on FCFs obtained from Morse potentials defined 
by spectroscopic constants available in the literature. Detailed population-dynamics 
studies are needed for a more rigorous selection. 

Appendix A.l. Diatomic molecular ions with transitions which excite a single electron 
from the HOMO to the LUMO 

Appendix A. 1.1. Group 3 hydrides and halides ScH"*" has an interesting ^A-^H system 
with very similar spectroscopic constants for both electronic states; however, the 
transition is deep into the IR (A = 6 /im) [7^. The ground state of YH"*" is a 
and the first excited state is a ^A. Spectroscopic constants re and Ue are similar but the 
transitions falls in the IR at ~ 3.3 /zm [75j. Finally, LaH"*" has a ^E+-^A system with 
similar spectroscopic constants and the transition is at ~ 5.1 /xm [76] . 
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Appendix A. 1.2. Group 13 hydrides and halides The open-shells from Group 13 and 
halogens are stable with ground states. Unfortunately, the ^11 potential curve is 
dissociative. This is common in the following examples: AlCl"*", AlF"*" [771 178] and 
GaCl"'". Although GaF"*" is believed to have strong ionic character [79] it also has a 
dissociative first exited state [80]. The hydrides BH"*" and AlH"*" seem to be suitable for 
direct laser cooling and are the subject of the present paper. Experimental spectroscopic 
data for GaH+ is not yet available [81] . 

Appendix A. 1.3. Bialkalis, alkali monohalides, alkali monohydrides and molecular 
ions formed by alkalis and Group 13 Li^ and LiNa"*" have and ^S"*" ground 

electronic states with only one valence electron, and they do not have valence correlation 
energy [82]. Their ground state dissociation energies are relatively high ( 10000 cm~^) 
and they do not seem to have dissociative states that couple non-adiabatically to the 
A^n-X^S"*" system. The lack of spectroscopic data on these systems stop us from making 
further predictions, but it would not be surprising if a few of these molecular ions have 
diagonal FCFs. These systems can not be ruled out as potential candidates. The same 
seems to be true for LiAl"*" with a bound ground state and a dissociation energy of 
roughly 1 eV [82]. LiB"*" has a weaker bond and the A and X states spectroscopic 
constants are different [83] . 

With respect to the alkali monohydrides, although they are technically open-shells 
with one unpaired electron spin, they have very weak bonds. This is expected as the 
same unpaired electron shares some density to participate in the bond (e.g., KH"*" [81] 
and NaH+ [85]) 

Appendix A. 1.4- Group 15 halides and hydrides PH+ and PF^ would require 
repumping of six and ten vibrational bands to scatter 2x10^ and 9x10^ photons, 
respectively. They do not have nuclear spin and their ground states are ^11. The first 
excited state of PH"*" is a and the vertical transition from the ground state is at 381 
nm. For the PF"*" case, the A^S-X^II optical transition is in the UV at 282 nm |86j . 

Appendix A. 2. Diatomic molecular ions with an unpaired electron in which an optical 
transition moves a hole to a bonding orbital 

Appendix A. 2.1. Group 17 hydrides Using Morse potentials derived from spectroscopic 
constants compiled by Radzigs and Smirnov [HS], we notice A^S-X^II transitions in 
molecules like HF"*", HCl"'" and HBr"*" to have very un-diagonal FCFs. However, a recent 
theoretical proposal to measure molecular parity violation [40j points out that HI"*" could 
have highly diagonal FCFs. This is attributed to the center of mass being very close 
to the iodine atom, resulting in the molecular orbit having mostly an atomic character. 
There is not enough spectroscopic information of excited states of HI+ to design a laser 
cooling experiment [87] . 
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Appendix A. 2. 2. Group 15 diatomic ions We use FCFs from spectroscopic constants 
compiled by Radzigs and Smirnov [86j to test N2' and P2'. The former could scatter 
11x10^ photons if irradiated by five lasers. These transitions fall in the near IR with 
the most energetic transition at 1091 nm. This homonuclear molecular ion has recently 
been loaded and sympathetically cooled in its rovibrational ground state [23j. Although 
P2' has a ^Uu ground state with similar re and Ue constants as its first excited state, 
the FCFs show that one would have to address more than 8 transitions with A > 4.6/im 
to scatter over lO'' photons. 

Appendix A. 2. 3. Group 5 and 6 halides These molecules have crowded excited state 

manifolds and the first electronic states high multiplicity. VF+ has a A^A — X^II system 

with very similar re and Ue constants and a vertical transition at ~ 7 /xm. If it is possible 

to cope with higher multiplicity (S=5), CrF+ might work as it has a A^Il — X^S"*" with 

an IR transition (A = 2 /xm) [88] . 
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